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Abstract 

A layer of uniaxial medium with large axial permittivity and per- 
meability can be used as a quarter-wave transformer with interesting 
properties. By increasing the transverse permittivity and permeabil- 
ity the transformer becomes a thin sheet. It is shown that the re- 
cently introduced SHDB boundary conditions, generalizing the soft- 
and-hard and DB conditions, realized by the interface of a skewon- 
axion medium, can be transformed to form a novel class of SHD'B' 
boundary conditions which generalizes the soft-and-hard and D'B' 
boundary conditions. Reflection of a plane wave from a planar SHD'B' 
boundary is considered by numerical examples revealing an interest- 
ing narrow beam with radical change of reflection for certain values of 
parameters and incidence angles. 

1 Introduction 



Various novelties in electromagnetic boundary conditions have been recently 
introduced and their physical and material realizations have been attempted. 
For example, the perfect electromagnetic conductor (PEMC) boundary con- 
ditions, 

nx(H + ME) = 0, n-(D-MB) = 0, (1) 
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involving one scalar parameter M were introduced in [1, 2] with material 
realizations suggested in [3, 4, 5]. Here, n denotes the unit vector normal 
to the boundary surface. As another example, the DB-boundary conditions 
requiring vanishing of the normal components of the D and B vectors were 
defined by [6, 7] 

n D = 0, n-B = 0. (2) 

They have found application in cloaking problems [8, 9, 10, 11, 12]. Some 
physical realizations for the DB conditions in terms of medium interfaces have 
been suggested in [13, 6, 14] and a more practical material realization was 
given in [15]. As a further example, boundary conditions requiring vanishing 
of normal derivatives of normal field components were defined in [16] by 

V • (nn • D) = 0, V • (nn ■ B) = 0, (3) 

which for the planar boundary z = are simplified to 

d z D z = D' Z = 0, d z B z = B' Z = 0. (4) 

Conditions (3) have been dubbed D'B' conditions. While a physical realiza- 
tion in terms of a medium interface has not yet been found, it was shown 
that the D'B' conditions can be obtained in terms of a layer of a certain 
medium upon the DB boundary [17, 18]. Generalizations of the DB and 
D'B' conditions have been discussed in [19, 20]. 

Finally, a novel set of conditions, generalizing those of the soft-and-hard 
boundary [21] and the DB boundary, was obtained in [22] when considering 
the interface of a certain class of skewon-axion media [23], also called IB 
media [24]. They were subsequently dubbed SHDB boundary conditions and 
expressed in the form 

An-B + u-E = 0, (5) 
AnD-uH = 0, (6) 

where u satisfying n ■ u = is a unit vector parallel to the boundary. For 
\A\ — > the conditions can be seen to reduce to those of the soft-and-hard 
boundary, u • E = 0, u • H = 0, while for \A\ — > oo they approach the 
DB conditions (2). Since all of the recently introduced boundary conditions 
(2), (3) and (5), (6) have been shown to be physically realizable in terms of 
medium interfaces, they cannot be considered just mathematical concepts. 
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However, there still is work to do to go one step further by realizing the 
medium parameters in terms of metamaterial structures. 

It is the purpose of the present paper to define a set of boundary con- 
ditions which is related to the SHDB boundary in the same way as the 
D'B' boundary is related to the DB boundary, i.e., a boundary which is a 
generalization of the soft-and-hard and D'B' boundaries. Because the D'B' 
boundary conditions could be produced by transforming the DB boundary 
conditions in terms of a layer of uniaxial medium [17], it appears reasonable 
to apply a similar idea for the SHDB boundary. 



2 Quarter-wave transformer 

Let us consider a uniaxial anisotropic medium defined by the medium equa- 
tions 

D = e t E t + e z u z E z , (7) 
B = ntHt + fJ, z u z H z , (8) 

where E t and H t are components transverse to the z axis. Assuming that 
the axial parameters become infinitely large, 

e z -> oo, fi z ->■ oo, (9) 

we must have E z — > and H z — > in the medium while D z and B z may have 
finite values. In practice it is sufficient that the parameters satisfy 

e 2 > IH, > IH- (10) 

Such a medium has been previously called by the name waveguiding medium 
[25], because a wave propagates along the axis of the medium like in a set of 
parallel waveguides. The axial components of the Maxwell equations are 

Uz-ViXE, = -juB z , (11) 
u z -V t xH t = juD z , (12) 

while the transverse components can be represented in the form 

d z E t = jk t r] t u z x Hi, (13) 
rj t d z B. t = -jk t u z x E t , (14) 
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with 

k t = uj^/J^Tt, r}t = vW e *- ( 15 ) 
Operating (13) and (14) by d z leads to the equations 

& + = ( °Y (16) 



Suppressing the dependence on x and y, the transverse fields must be of the 
form 

E t (z) \ ( K, \ jktz f K, » jh . 



h,w y - v h, + j e + { Hr ) e ' (17) 

Actually, (E+,H+) and (E^,H^) may have any dependence on x and y. 
Inserting (17) in (13) and (14) and equating components with similar z de- 
pendence, we obtain the following relations for the field amplitude vectors: 

E+ = -u z x Vt n+ (18) 
E- = u 2 x^H-. (19) 

Let us now consider the fields in the waveguiding medium at z = as 
related to fields at z = — d when choosing 

d = ^ = j, exp(±jk t d) = ±j. (20) 

The relations can be expressed as 

Et(0) = jr] t u z xU t (-d), (21) 
7^H t (0) = -ju 2 xE t (-4 (22) 

The corresponding relations for the axial field components are obtained by 
substituting (17) in (11) and (12), and the same differentiated by d z . Skipping 
the details, the results can be represented as 

DM \ 1 ( D' z (-d) \ 

D ' M ) ~ -k t ( Dz{ - d) \ (24) 
B'M) ~ kt \B z (-d)J- W 

Equations (21), (22), (23) and (24) represent relations between fields at 
the two planes in the waveguiding medium. A layer of quarter wavelength of 
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such a medium transforms boundary conditions at z — —d to other boundary 
conditions at z — 0. By letting the transverse parameters grow in magnitude 
with respect to the parameters e,/x of the exterior medium, e t /e — > oo and 
lit I '/x oo, the layer becomes a thin sheet. In such a case there exist two 
different orders of large parameters, e z e t 3> e, \i t ^> /x t ^ I 1 - 



-d 





/ e 
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Figure 1: Quarter- wave transformer is a layer of uniaxial medium with large 
axial parameters, \i z ^> fi t and e 2 ^> e t . It can be used for transforming 
boundary conditions at z — —d to other boundary conditions at z = 0. For 
[At ^> n and e t ^> e the layer becomes a thin sheet. 



One must note that, while E t , H t , D z and i? x are continuous through the 
interfaces, D' z and B' z are not. Thus, (23) and (24) are valid only when D' z 
and B' z are taken at the interior side of the interface. Denoting the interior 
side of the interfaces by Qi and the exterior side by () e , the relations can be 
obtained from 



V-B, = 8 z B zi + fxtVt ■ Hti = 0, (25) 

V'Dj = d z D zl + e t V t ■ E tl = 0, (26) 

V B e = d z B ze + fiV t ■ H ie = 0, (27) 

V D e = d z D ze + eV t -E te = 0. (28) 



Let us assume that the half space z > is isotropic with parameters e, /x 
while the other half space z < —d is taken care of by boundary conditions 
at z = —d e . Due to continuity of Vt • E t and Vf • H t through the interface 
z = 0, we obtain the relations 

fiB'M) = fi t B' z (0 e ), eD'M) = e t D' z (0 e ). (29) 



Assuming for simplicity /x t //x = e t /e = fct/fc, with k = uj^/JIe, the conditions 
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(23) and (24) can be written as 

D z (O e ) \ 1 / D'^-di) 



B z (O e ) I ~ kt \ B>(-d t ) ! (30) 



B'M) - M B z (-d t ) 



(31) 



3 SHD'B' boundary conditions 

The case of obtaining D'B' conditions (4) by transforming the DB conditions 
is now obvious from (31). Also, other cases considered in [17] are directly 
obtained. Let us now consider the SHDB conditions (5), (6) at z — —d e , by 
assuming u = u^,, 

AB z (-d e ) + E x (-d e ) = (32) 
AD z (-d e ) - H x (-d e ) = 0. (33) 

Since these field quantities are continuous through the interface z = —d, 
we can replace d e by di and apply (21), (22) and (31). Thus, the SHDB 
conditions are transformed through the quarter-wave transformer to the form 

AB' g (0 e ) - jufjiH y (0 e ) = 0, (34) 
AD' z (0 e ) — jueE y (0 e ) = 0. (35) 

Because the boundary conditions (34) and (35) correspond to a certain physi- 
cal setup, a layer of waveguiding medium upon a half space of certain medium 
yielding the SHDB conditions, they have a firm physical significance. It is 
proposed that they will be called SHD'B' boundary conditions. From the 
analogy of conditions (3) and (5), (6), the more general form for the SHD'B' 
boundary conditions is 

AV ■ (nn • B) + v ■ (V x E) = 0, (36) 
AV ■ (nn • D) - v • (V x H) = 0, (37) 

where v = n x u for u defined by (5), (6). In this form the SHD'B' conditions 
are valid for more generally than for planar boundaries. Making the trans- 
former a thin sheet, and because its transforming property is local, it can 
be shaped to follow any smooth surface (with radius of curvature assumed 
much larger than the thickness of the sheet). 
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The two extreme values of the parameter A yield two known special cases: 
\A\ — > corresponds to the conditions of the soft-and hard boundary while 
\A\ — > oo corresponds to those of the D'B' boundary. Since in the form (36), 
(37) the boundary conditions do not involve the frequency, they are valid to 
any time dependence like the SHDB conditions (5), (6). Furthermore, both 
sets of boundary conditions are self dual. In fact, making the simple duality 
transformation [26] 

E d = j V H, U d = E/ 3V , (38) 

with rj = y/Ji/e, the pair of conditions (36), (37) can be seen to remain 
invariant. Thus, just like the SHDB boundary conditions (5), (6), the SHD'B' 
boundary conditions are self dual. 

4 Plane- wave reflection from SHD'B' plane 

To have an idea of the properties of the SHD'B' boundary, let us consider 
the basic problem of plane-wave reflection from a SHD'B' plane z — 0. The 
incident plane wave in the isotropic half space z > is defined by 

E*(r) = E i exp(-jk i • r), (39) 
H*(r) = H l exp(-jV • r), (40) 

and the field reflected from the SHD'B' is denoted by 

E r (r) = E r exp(-jk r • r), (41) 
H r (r) = H r exp(-jk r • r). (42) 



The different k vectors are 



k l = -k z u z + k t , 
k r = k z u z + k t , 

k( — k x ii x + kylly. 



(43) 
(44) 
(45) 



Applying the Maxwell equations for the fields of the incident and reflected 
plane waves, 



kri& = k* x E\ kr]B. r = k r x E r , (46) 

kW = -k* x r?H\ kW = -k r x r/H r , (47) 
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and the orthogonality conditions 

k'-E* = k r -E r = 0, (48) 
k'-H* = k r -H r = 0, (49) 

relations between the tangential field components can be compactly expressed 
as 

77Hj = -J.Ej, ^=J.E?, (50) 

Ei = J.^Hi, E r t = -J. V W t . (51) 

The 2D dyadic 1 defined by 

= J = -^rUzX(k% + k t k t ), (52) 
resembles the imaginary unit because of its properties [6] 

I 2 = -I, T 1 = -1, (53) 

which can be easily verified. 

The SHD'B' conditions (34), (35) can be expressed in terms of a dimen- 
sionless parameter C defined by 

Ah 

C= — = A^/JTe, (54) 
ui 

in a form involving the tangential field components as 

Cd x H x (p) + (Cdy + jh)H v (0) = 0, (55) 
Cd x E x (0) + (Cd y + jk)E y (0) = 0. (56) 

For the plane wave fields, (55), (56) become 

Ck x H x (0) + (Ck y - k)H v (0) = 0, (57) 

Ck x E x (0) + (Dk y — k)E y (0) = 0. (58) 

Defining the vector 

c = Ck x u x + (Cky — k)u y = Ck t — ku y , (59) 
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parallel to the boundary, the conditions (57), (58) can be expressed for the 
plane-wave fields in the simple form 

c-H(O) = 0, (60) 
c-E(O) = 0. (61) 

These conditions have an appearance similar to that of the soft-and-hard 
conditions [21]. However, here the vector c depends on the k* vector of the 
incident wave. 

In terms of the amplitude vectors of the incident and reflected field com- 
ponents, the SHD'B' conditions (60), (61) have the form 

c(Hj + H[) = 0, (62) 
c(Ei + E[) = 0. (63) 

Applying (50), (51), the conditions (62) and (63) can be respectively rewrit- 
ten as 

b.(Ej-EJ) = 0, (64) 
b.(Hj-HJ) = 0, (65) 

where b is another vector parallel to the boundary, 

b = c • J = -^—(k 2 z c x u z - (u 2 • c x k t )k t ). (66) 



kk 



z 



Conversely, we can write 

c = -b • I (67) 

In the general case the vectors b and c are linearly independent, i.e., they 
satisfy 

bxc-u^-j-fcxk'jVO, (68) 

whence they make a vector basis in the boundary plane. The vectors b', c' 
defined by 

V = (69) 

C = (70) 
u, • b x c 
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form the reciprocal basis satisfying 

b ■ c ' = c • b' = 0, b • b' = c • c' = 1. (71) 
The 2D unit dyadic can be expressed as 

\ = b'b + c'c. (72) 
From (66) and (67) we can write the expansion 

J = \ . J = c 'b - b'c, (73) 
whence the reciprocal basis vectors satisfy the relations 

c' = J • b', b' = -l- c'. (74) 

5 Eigenfields 

The 2D reflection dyadic R can be defined by the relation between the tan- 
gential field components as 

E? = I • Ej. (75) 

Applying (63), (64) and (72), the reflection dyadic in (75) can be constructed 

as 

W = b'b -c'c. (76) 
Uj, x (cb + be) 



u z -cxk t 
Defining the angles of incidence 6, ip by 



(77) 



k x = /csin^cosv?, (78) 

k y = A; sin 6* sin (79) 

k z = kcos6, (80) 

the components of the expansion 

R \1 X U X R XX ~\- \l X 1\yR X y -\- WyW X Ry X ~\~ WyWy Ryy ^ ( 8 1 ) 
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can be expressed as 

R xx = -^-((cos 2 9 + sin 2 9 cos 2 ip)(l — 2Csin#sin ip) 

-C 2 sin 2 9 cos 2 9 cos 2ip) = -R yy , (82) 

R xy = ^■(sin 2 ^sin2v?(l - C 2 cos 2 #) 

+2C sin 9 cos y?(cos 2 9 - sin 2 9 sin 2 ip)), (83) 

R yx = -^-(2(7 sin # cos y9 (cos 2 # + sin 2 # cos 2 y?) 

-<7 2 sin 2 #cos 2 #sin2y9), (84) 
A = cos 2 ip + (Csin6» - siny?) 2 cos 2 6>. (85) 

One can easily show that the coefficients R xx , R yy are symmetric, and R xy , R yx 
are antisymmetric, in the variable n/2 — (p. 

Since the eigenproblem for the reflection dyadic (76) has the solutions 

W • b' = b', R • c' = -c', (86) 

the eigenfields corresponding to the eigenvalues +1 and —1 are respectively 

Ei + = £ + b' = E[ + , (87) 
E\_ = E^c' = -W t _. (88) 

The corresponding magnetic eigenfields are obtained by applying (50), 

V W t+ = -l-El + = -E + c' = - V W t+ , (89) 
V W t _ = -3-El = E_b' = riW t _. (90) 

Since the total eigenfields satisfy 

(EL + EI), = 0, (91) 

(hv + h;), = o, (92) 

the — field is reflected as from a PEC plane and, the + field, as from a 
PMC plane. A corresponding property was previously found for the SHDB 
boundary [22]. 
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Applying again k l • E!j_ = and k* • H± = 0, expressions for the total 
eigenflelds can be reduced to 





= E + 


E!_ 


= 




= E + 




= E_ 



cxk 1 
k« • b x c ' 

k* x b 
k« • b x c ' 

bxk ! 
k* • b x c ' 

cxk 1 



(93) 
(94) 
(95) 
(96) 



k'-bxc 

The field amplitudes are obtained from the total field as 

E + = b ■ E\ E- = c ■ E\ (97) 

6 Special cases 

Let us consider some special cases of the SHD'B' medium and the incident 
plane wave. 

Soft-and-hard boundary 

For the case C = the SHD'B' boundary reduces to the soft-and-hard bound- 
ary with boundary conditions 

u y • E = 0, • H = 0. (98) 

Because (59) yields c = — ku y , and from (93), (96) both E+ and H*_ have 
the polarization Uj, x k ! , they satisfy 

u y • E\ = 0, u y - Hi = 0. (99) 

Thus, the + and — waves can be respectively called TE^ and TlVLj waves. 
The reflection dyadic is simplified to 

= sin 2 6 sin 2(p , , 

R = u x u, - + — 2~~p. ^— j^u z u r (100) 

sin U cos^ (f + cos^ v 
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D'B' boundary 

For C — > oo we have c — > Ch t , whence both and H!_ have the polarization 
k t x k ! = k z k t x u 2 . In this case the eigenfields satisfy 

u 2 ■ E\ = 0, u 2 ■ IV = 0. (101) 

The + and — waves can now be respectively called TE 2 and TM 2 waves. In 
this and the previous case field decomposition in TE and TM components is 
independent of the k* vector, whence it is valid for any fields. The reflection 
dyadic can be written in this case as 

R = - cos 2y?(u ;c u ;E - u y u y ) - sin 2(f(u x u y + u y u x ). (102) 

Special angles of incidence 

Let us consider special cases of the vector k\ 



In the case u x ■ k* = k x = 0, or <p = n/2 the reflection dyadic (76) has 
the simple form 

R = u,u x - UyUy. (103) 

For a wave with this incidence the boundary appears as a hard surface 
[21]. 

The same result is obtained for the case k t = 0, or 9 = 0. 

In the case • k l = k y = 0, or (p = 0, the reflection dyadic takes the 
form, 

1 - C^sin^cos 2 ^ 
1 + C sin 9 cos 9 



2C sin 6>(cos 2 9u x u y + u y u x ) 
l + C 2 sin 2 flcos 2 # 



This expression coincides with the ones above for C = and C — > oo. 
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7 Numerical examples 



Let us consider as representative examples the reflection coefficients i?y for 
some incident plane waves and some values of the SHD'B' boundary param- 
eter C. Figure 2 illustrates the four reflections coefficients as functions of 
the incidence angle 9 when the plane of incidence is the xz plane (ip = 0). 
The figure shows clearly that R xx = —R yy for all parameter values of 9 and 
C. Another interesting observation is that for C = 2 and 9 = 45°, the 
co-polarized reflection vanishes. In other words the surface acts as a twist 
polarizer. 



fix 



\ 



/ J 



30 60 90 




30 60 90 



6 




30 60 90 



R. 



-1 



SO 60\ 90 

/ \ i 



Figure 2: The four components of the reflection matrix as function of the 
incidence angle 9 with three values for the parameter C: solid line (C = 0.5), 
long-dashed line (C = 2), short-dashed line (C = 20). The azimuth angle is 
tp = 0. 



Due to the two-dimensional anisotropy of the surface, the reflection char- 
acteristics should depends on the azimuthal angle ip. This is displayed in 
Figure 3 where the co-polarized component R xx ((p) and the cross-polarized 
component R xy ((p) are depicted for 9 = 0.3 and C = 3. The curves show 
that R xx is symmetric with respect to the incidence ±y-axis (<p = 90° and 
ip = 270°) whereas i? xy is antisymmetric. Similarly, R yy turns out to be 
symmetric and i? yx antisymmetric. The effect of the narrow-beamed radical 
change in reflection occurs when both cos ip and C sin 9 — 1 have small values. 
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In such a case the reflection coefficient can be approximated by 

_ (Csinfl- l) 2 cos 2 g - cosV 
(C sin # — 1J Z cos^ y + cos z ip 

At cos<y9 = we have R xx = 1 (exactly) while the nulls of R xx appear at 
cos <p ~ ±(C sin 9 — 1) cos 9. 




Figure 3: The components R^ (solid line) and i? xy (dashed line) of the 
reflection matrix (75) as function of the azimuth angle ip for incidence angle 
9 = 0.3 and the parameter C = 3. 

Figure 4 shows how the reflection approaches that of a hard surface when 
the incidence plane comes closer to yz plane, if — > tt/2. For the hard surface 
has i? xx = +1 independently of the incidence angle, as can be seen from 
the figure. Again, an interesting narrow-beamed change of reflection at 9 = 
arcsin(l/C) ~ 0.524 = 30° can be observed. 



8 Conclusion 

A novel set of electromagnetic boundary conditions was defined in the present 
paper. Because the soft-and-hard conditions and the D'B' conditions appear 
as special cases of the boundary conditions, the generalization was dubbed 
SHD'B' conditions. The present paper is a continuation to a previous one 
[22] generalizing the soft-and-hard and DB conditions to what were called 
the SHDB conditions. Since the D'B' conditions can be realized in terms of 
a medium layer transforming the DB conditions, the same idea was applied 
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Figure 4: R xx as function of the incidence angle 6 when the plane of incidence 
is close to yz plane (ip = vr/2). The azimuth angle is <p = n/2 — 0.01 (solid 
line), if = 7r/2 — 0.1 (long-dashed line), and (p = n/2 — 0.2 (short-dashed 
line). The boundary parameter has the value C — 2. 

to the realization of the SHD'B' conditions. Making the layer a thin sheet by 
proper choice of the parameters makes it possible to apply the transformation 
for the SHDB boundary of any form. To see basic properties of the boundary, 
reflection of a plane wave of a planar SHD'B' boundary was considered. An- 
alytic expressions for the reflection dyadic and its eigenfields were derived. It 
was shown that for the two eigenpolarizations the SHD'B' boundary appears 
as the PEC or the PMC boundary with the respective reflection coefficients 
+1 and —1. Numerical computations have shown that, for certain values of 
the parameter C and angle of incidence, the field reflected from the SHD'B' 
boundary is totally cross polarized. Also, there exists an interesting narrow- 
beam region for the incident wave where the reflection makes a radical change 
in magnitude. Such a propery could have an application similar to the corner 
reflector in defining a given direction in space. This may serve as a reason to 
study the properties further and attempt to find a metamaterial realization 
for the SHD'B' boundary. 



References 



[1] Lindell, I. V. and A. Sihvola, "Perfect electromagnetic conductor", J. 
Electro. Waves Appl. Vol. 19, No. 7, 861-869, 2005. 



17 



[2] Sihvola, A. and I.V. Lindell, "Perfect electromagnetic conductor 
medium", Ann. Phys. (Berlin), Vol. 17, No. 9-10, pp.787-802, 2008. 

[3] Lindell, I. V. and A. Sihvola, "Realization of the PEMC boundary," 
IEEE Trans. Antennas Propag., vol.53, no. 9, pp. 3012-3018, September 
2005. 

[4] Shahvarpour, A., T. Kodera, A. Parsa, and C. Caloz, "Arbitrary elec- 
tromagnetic conductor boundaries using Faraday rotation in a grounded 
ferrite slab," IEEE Trans. Microwave Theory Tech., vol.58, no.ll, 
pp.2781-2793, 2010. 

[5] El-Maghrabi, H.M., A.M. Attiya and E.A. Hashish, "Design of a per- 
fect electromagnetic conductor (PEMC) boundary by using periodic 
patches," PIER M, vol.16, pp. 159-169, 2011. 

[6] Lindell, I. V. and A.H. Sihvola, "Uniaxial IB-medium interface and novel 
boundary conditions," IEEE Trans. Antennas Propagat., Vol. 57, No. 3, 
pp.694-700, March 2009. 

[7] Lindell, I. V. and A. Sihvola: "Electromagnetic boundary condition and 
its realization with anisotropic metamaterial," Phys. Rev. E, Vol. 79, 
No. 2, 026604 (7 pages), 2009. 

[8] Zhang, B, H. Chen, B.-I. Wu and J. A. Kong, "Extraordinary surface 
voltage effect in the invisibility cloak with an active device inside," Phys. 
Rev. Lett, Vol. 100, 063904 (4 pages), February 15, 2008. 

[9] Yaghjian. A. D. and S. Maci, "Alternative derivation of electromagnetic 
cloaks and concentrators," New J. Phys., Vol. 10, 115022 (29 pages), 
2008. Corrigendum, ibid, Vol. 11, 039802 (1 page), 2009. 

[10] Weder, R. "The boundary conditions for point transformed electromag- 
netic invisible cloaks," J. Phys. A, Vol. 41, 415401 (17 pages), 2008. 

[11] Kildal, P.-S., "Fundamental properties of canonical soft and hard sur- 
faces, perfect magnetic conductors and the newly introduced DB surface 
and their relation to different practical applications including cloaking," 
ICEAA'09, Turin, Italy, August 2009, pp.607-610. 



18 



[12] Martini, E., S. Maci and A.D. Yaghjian, "DB boundary conditions at 
the inner surface of an arbitrarily shaped cloak," EuCAP 2011, Rome, 
Italy, April 2011, pp.3609-3611. 

[13] Rumsey, V. R. "Some new forms of Huygens' principle," IRE Trans. 
Antennas Propagat., Vol. 7, Special supplement, pp.S103-S116, 1959. 

[14] Lindell, I.V. and A. Sihvola, "Simple skewon medium realization of DB 
boundary conditions", PIER Letters, vol.30, pp.29-39, 2012. 

[15] Zaluski, D., D. Muha and S. Hrabar, "DB boundary based on reso- 
nant metamaterial inclusions," Metamaterials '201 1 , Barcelona, October 
2011, pp.820-822. 

[16] Lindell, I. V., and A. Sihvola, "Electromagnetic boundary conditions 
defined in terms of normal field components," Trans. IEEE Antennas 
Propag., Vol. 58, No. 4, pp. 1128-1 135, April 2010. 

[17] Lindell, I.V., A. Sihvola, L. Bergamin and A. Favaro: Realization of 
the D'B' boundary condition, IEEE Antennas Wireless Propag. Lett., 
vol.10, pp.643-646, 2011. 

[18] Lindell, I.V., J. Markkanen, A. Sihvola and P. Yla-Oijala, "Realization 
of spherical D'B' boundary by a layer of wave-guiding medium;" Meta- 
materials, vol.5, pp.149-154, 2011. 

[19] Lindell, I.V., H. Wallen and A. Sihvola, "General electromagnetic 
boundary conditions involving normal field components," IEEE Anten- 
nas Wireless Propag. Lett., vol.8, pp. 877-880, 2009. 

[20] Wallen, H., I.V. Lindell and A. Sihvola, "Mixed- impedance boundary 
conditions," IEEE Trans. Antennas Propag., vol.59, no. 5, pp. 1580-1586, 
May 2011. 

[21] Kildal, P.-S. "Definition of artificially soft and hard surfaces for electro- 
magnetic waves," Electron. Lett., vol.24, pp. 168-170, 1988. 

[22] Lindell, I.V., and A. Sihvola, " Skewon- axion medium and soft-and- 
hard/DB boundary conditions," IEEE Trans. Antennas Propag., sub- 
mitted. 



19 



[23] Hehl, F. W. and Yu. Obukhov, Foundations of Classical Electrodynam- 
ics, Birkhauser, Boston, 2004. 

[24] Lindell, I. V. "The class of bi-anisotropic IB media," Prog. Electromag. 
Res., Vol. 57, pp.1-18, 2006. 

[25] Lindell, I.V. and A.H. Sihvola: "Realization of impedance boundary," 
IEEE Trans. Antennas Propag., vol.54, no. 12, December 2006, pp. 3669- 
3676. 



[26] Lindell, I. V., Methods for Electromagnetic Field Analysis, Oxford Uni- 
versity Press, 1992. 



